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Abstract 

The following principle of minimum energy may be a powerful substitute to the dynamical 
perturbation method, when the latter is hard to apply. Fluid elements of self-gravitating 
barotropic flows, whose vortex lines extend to the boundary of the fluid, are labelled in 
such a way that any change of trial configurations automatically preserves mass and circu- 
lation. The velocity field is given by a mass conserving Clebsch representation. With three 
independent Lagrangian functions, the total energy is stationary for all small variations 
about a flow with fixed linear and angular momenta provided Euler's equations for steady 
motion are satisfied. Thus, steady fiows are stable if their energy is minimum. Since en- 
ergy is here minimized subject to having local and global contants of the motion fixed, 
stability limits obtained that way are expected to be close to limits given by dynamical 
perturbation methods. Moreover, the stability limits are with respect to arbitrary, not 
necessary small, perturbations. A weaker form of the energy principle is also given which 
may be easier to apply. 

The Lagrangian functional, with the same three Lagrange variables is stationary for 
the fully time dependent Euler equations. It follows that the principle of minimum energy 
gives stability conditions that are both necessary and sufficient if terms linear in time 
derivatives (gyroscopic terms) are absent from the Lagrangian. The gyroscopic term for 
small deviations around steady fiows is given explicitly. 

Key words: Energy variational principle; Self-gravitating systems; Stability of fluids. 
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1. Introduction 

The method of small perturbations and analyses of eigenmodes are difficult to apply- 
to astrophysical models that are not at least axisymmetric (Binney & Tremaine 1987). 
Moreover, it gives good indications on how instabilities set in, but once a configuration 
has become unstable, the method is at a loss to say whether it tends towards a new steady 
figure or will oscillate like a pendulum or become chaotic. In this era of supercomputers, 
one might think that a dynamical follow-up is routine work, but this is not the case, 
yet. Moreover, as a rule, analytical methods give more insight than purely numerical ones. 
Therefore, there are still reasons for perfecting other methods like the principle of minimum 
energy which played important roles in stability theory [see in particular Antonov (1962), 
Lebovitz (1965) and Chandrasekhar's (1969) works]. 

We shall be concerned with ideal fluids and, for definiteness, we limit ourselves to 
barotropic fluids. In our Lagrange variables, non-barotropic fluids are more or less a 
particular case than a different one and may be treated along the same lines. We shall not 
do this here. 

From the pioneering work of Arnold (1966), see also Lynden-Bell & Katz (1981) on 
compressible flows (the paper will be referred to below as LBK81) we know that among all 
steady flows that are mass preserving and isocirculational, those whose energy is stationary 
for small variations satisfy Euler's equations. By mass preserving we mean that in an 
arbitrary virtual displacement of a small element of fluid there is no loss of mass and by 
isocirculational we mean that the flux of the vortex lines, in any direction through a small 
area, is the same for the displaced element of fluid. 
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It follows that if a steady configuration has minimum energy, it is stable with respect 
to perturbations that are not necessarily small. In a perturbation analysis, perturbations 
are always assumed to be small. Thus, while there may be stable configurations with 
respect to small perturbations for which the energy is not a minimum, the stability limits 
given by the energy principle are stronger than those of a perturbation analysis. 

This principle of minimum energy can be made even stronger by imposing further, 
as constraints, global invariants of the motions. In that way, virtual displacements may 
eventually be restricted as much as real motions. If, by varying some parameter, a series 
of stable configurations ceases to have minimum energy with such severe constraints, the 
limit may be quite close to the limit one would find in a perturbation analysis. 

Another way of using the "maximally constrained" energy consists in using a computer 
to find a new minimum when the energy principle fails to indicate stability. Since time is 
not involved, calculations have one dimension less than the dynamical equations ; this may 
be of great computational help. The new configuration, if any, will be one in which the 
unstable model is likely to settle. The absence of a new steady configuration may indicate 
that no new real stable configuration exists. 

A key point in applying this energy principle is to make an appropriate parametrization 
of the fluid elements; this includes the topology of the flows that must be conserved as 
well. The parametrization used here has been found in LBK81 where labels are attached to 
vortex lines rather than to fluid elements. Such a labelling leads naturally to a particular 
Clebsch representation for the velocity of the flow. A Lagrangian and Energy Principle 
in terms of three independent functions was developped in Katz and Lynden-Bell(1985) 
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(For a formal generalization of that approach, see Simo, Lewis and Marsden 1991). The 
method contains no Lagrange multiplier as is usually done [Serrin 1959, Lin 1963, Seliger 
and Whitham 1968] as such multipliers are not very usuful in second variations. 

The new elements in this work, as compared to LBK 81 and KLB 85, are as follows. We 
specify completely the labelling in a class of flows. We use a new set of Lagrangian functions 
to prove energy principles for steady motions. We show the uniqueness of our choice and 
fix whatever freedom remains in positioning the coordinates. The principle of minimum 
energy is completed with the proof that our variables are true Lagrangian variables. That 
is, Euler's dynamical equations are obtained by varying the Action. If the Lagrangian's 
kinetic energy is purely quadratic in time derivatives (no gyroscopic term), the principle 
of energy minimum is a necessary and suflftcient condition for stability. The form of the 
terms linear in time derivatives gives thus important information and have therefore been 
worked out explicitly. The formulation follows closely classical mechanics. Applications 
to two dimensional flows and in particular to MacLaurin disks, whose dynamical stability 
limits are well known, give a perfect illustration of the power of our method. To maintain 
this paper within reasonable limits, we delay the application to an accompanying paper IL 

2. Stationary Barotropic Self-Gravitating Flows 

It will be useful to write the equations of motion first to show our essentially standard 
notations: f — {x^) = {x,y,z) [k,l,m,n — 1,2,3] for the positions, W for the steady 
velocity field in inertial coordinates, p the density, P{p) the pressure, h{p) = J dP/p the 
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specific enthalpy and $(< 0) the gravitational potential, 

^ = -G I -^^Q-d^x'. (2.1) 
J \r — r I 

Euler's equations for flows that are steady in coordinates with uniform velocity b and 
angular velocity* velocity U may be written 

O = (U ■V)W + ilcXW + Vih + ^) = 0, (2.2) 

where, by definition, 

W = U + b + ncXf=U + f]c. (2.3) 
Another useful form of Euler's equation is 

d = a;xt/ + VA = 0, (2.4) 

where 

a; = rot (2.5) 

and 

A=^(t/2-77^) + /i + $. (2.6) 
With (2.4) we also have the equation of mass conservation 

U = V- {pU) = 0. (2.7) 

The equation of circulation conservation follows from (2.4): 

rot (ujxU) = 0. (2.8) 



* The index c of Qc referes to rotating coordinates and is just here to distinguish it from 

— * 

Q with no index, commonly used for angular velocities in galactic disks or rotating ste 
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The boundary conditions are those of a self-gravitating flow in free space, the density and 
the pressure go to zero and the velocity of sound goes to zero as well: 



P|. = ^ 

dp 



= 0. (2.9) 

p=0 



Notice that (2.7) with (2.9) implies 

^•VpU = 0, (2.10) 

— * 

i.e., U is in the tangent plane of the surface of the fluid. Apart from satisfying these 
boundary conditions, physical quantities are assumed to be bounded everywhere. 

In time dependent flows, in addition to mass conservation and conservation of circu- 
lation, the total mass, the linear and the angular momentum are conserved; the center of 
mass moves with a uniform velocity or is at rest. The conserved global quantities are 

M = j pd^x, P = Jwpd^x, J = JrxWpd^x. (2.11) 

It is always possible and worthwhile to take 

P = 0. (2.12) 

What this has to do with steady flows will soon become clear. 
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3. Mass Preserving and Isocirculational Labelling 

3.1. Labelling of Fluid Elements When the Vortex Lines Extend to the Boundary 

Conservation of vorticity implies conservation of the topology of the vortex lines. 
Flows with given vorticity have thus also given topologies of tJ-lines. For definiteness we 
shall consider flows with the same topology as in LBK81 (see figure 1). There, vortex lines 
were labelled by three parameters: the load A, the metage /j, and an "angle" p. 

Surfaces of constant load A are surfaces of constant mass per unit vortex strength in 
a narrow tube of vortex lines: if dM is the mass in a tube with vorticity flux dC, 

A=^=r^^'. (3.1) 



^Bottom ^ 

The integration is along a;-lines. The surfaces A = constant are embedded "cylinders" as 
shown in figure 1. The cylinders may be parametrized in any way we want, say a — a (A) = 
const., but the particular parametrization in which a (A) is proportional to the circulation 
on a closed contour around A = const, has very special property (see below). 

The metage /x, is defined as another family of surfaces with the same integral as in 
(3.1) taken up to some "red mark" on the tD-line: 

/■Red mark 

fM= f^dl. (3.2) 

isottom ^ 

Surfaces of constant ji cut accross A-surfaces. If the red mark is on the bottom, = 0. If 
it is at the top, fj, = X. 

The angular variable /? is defined by a family of surfaces of tD-lines as well ; these are 
"cuts" hanging on a "central line" (see figure 1). Thus, by definition, 

cl; • VcK = 0, a; • V/? = and u-V^ = p. (3.3) 



The parametrization of (3 may be so chosen that 

J} = Va X V/3. (3.4) 

This defines (5 for given uj and a up to an additive single valued function of a, B{a). This 

B{a) is associated with the freedom to take any cut as /? = 0. Having chosen a,P,fj,, It 
also follows from (3.3) that we must have 

p = Va X V/3 • V// (3.5) 

or that 

d{x, y, z) 



p=|^l^. (3.6) 



and thus 

5(tt,/3,/u) 



^ (3.7) 



a(x,2/,2;) 

everywhere except on the surface of the fluid since p\s = 0, 

p\s=Cj-Vi^\s = 0. (3.8) 

With (3.4) we see that the velocity field has now a Clebsch form 

W = aVp + Vi^. (3.9) 

The function u we define by the condition that mass be preserved i.e. by equation (2.7). 
The function u is thus defined by an elliptic equation since (2.7) can be written 

V- (pW) = V- [p(-aV/? + 77c)]. (3.10) 

Since p\s = 0, a regular single valued solution z/(r) (up to a constant) must be unique when 
it exists. 



The main tool in practical calculations will thereby be the Green function G{f,f'), 
solution of 

V- (pVG) =5^(f-f '). (3.11) 



3.2. Uniqueness of the Labelling 

Consider now the effect of a reparametrization of the fluid elements in a flow with 
given p{r), W{r) and thus uj{f) and given X{a). Let a,P,p, be another labelling. Surfaces 
of constant load may be reparametrized in any way: a — > a{a). Since 

—* —* ~ dd —* — * ~ 
= Va X V/3 = Va X V/3 = — Va x V(3 (3.12) 

da 



We must change (3 to 



P=^(3 + B{a). (3.13) 
da 



With both /i = /i = on the bottom of the fluid, the parametrization of fj, for given p and 
uj is uniquely defined 

A = /X. (3.14) 
The corresponding change in u is accordingly defined by 

W = aVP + Wu = aWP + WD (3.15) 

or 

V{u-u)^aV(3-aV(3. (3.16) 

From (3.13) it follows that 

i> = iy- (a-^ - a)l3 - / adB{a) + Const. (3.17) 
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a is single valued but j3 is not. If, however, we define and keep v single valued, any 
reparametrization that preserves the single valuedness of v must satisfy the condition 

da 

OL— — a = or 6l = la^ (3.18) 
da 

where I is a constant. The parametrization for which v is single valued is easily found. On 
the one hand C(A) is given by the flux of uj through a surface with boundary a = constant. 

C(A) = j Jo-dS = J Jdadp^ j [(3]da, (3.19) 

where [f3] is the value of the discontinuity of (3 and etc is the value of a on the central line. 
On the other hand C(A) is also equal to the circulation of W along a contour on a a = 
constant surface: 

C(A) = ([ W-df= ([ adp = a[P]. (3.20) 

J a J a. 



So, (3.19) and (3.20) give 



§^M^„l.„f (3.., 



which means that 



[P] = q and a = ^C(A), q = canst. (3.22) 

The parametrization of a by the circulation along A-tubes is unique except again, as in 
(3.18), for the multiplication constant g; with q — 27r, the domains of the Lagrange variables 
a, P, fi is thus 

< a < CKM = CMAxf^TT, B{a) <p< B{a) + 27r, < < \{a), (3.23) 

where a — -^^^ must he a given function of the flow. C{\) is similar to the amplitude 
J of the angular momentum. The only remaining arbitrariness of the parametrization is 

11 



B{a) associated with the cut where /? = 0. We have thus found that the parametrization 
a = C/2iT of LBK81 is the only one (up to q) that insures single valued z^'s in a Clebsch 
representation. 

3. 3. The Fixation of Coordinates and of the Cut (3 — 

Consider the trial configuration of figure 1. We are free to position the axis of coor- 
dinates in the simplest way; we shall make the following choice. The central line has two 
points fj, = 0, fj, = Xc', we use this as the z axis. The bottom of the central line is the origin 
of the coordinates. Unless the configuration is axially symmetrical, we have several ways, 
for orienting the xy coodinates. For instance if the yz plane cuts the A = line at P, 
the orientation may be chosen so that yp be maximum. There may be several extrema of 
y^ and it does not matter which one we take. We shall only compare configurations with 
small differences. The cut /3 = may be the ruled surface generated by lines parallel to 
the X axis and sliding on the central string. 

In this way, we fix uniquely the relative positions of trial configurations with respect 
to the coordinates and the parametrization of their points is: 



3.4- Proof That the Labelling Is Mass and Circulation Preserving 

Consider a fiuid element labeled {a, P, jj) = {a^) with the coordinates a;, j/, z or r. Any 
displacement from r to r + Ar of that fiuid element (Act = A/3 = A// = 0) is obviously 
mass preserving since according to (3.5) 



0<a<aM, 0</3<27r, < < A(a) 



(3.24) 



/^{pd^x) = Ad^a = 



(3.25) 
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This shows incidentally that if we set 

Ar = i{r), then Ap = -pV-C (3.26) 

Similarly from (3.4), the flux of the vorticity through any surface element dS parametrized 
by a, (3 is 

Co -(18 = dad(3 (3.27) 

which shows again that 

A(a; ■dS)^Q (3.28) 

or (see Arnold 1966, LBK1981) 

AcD + rot(cD X f) = {(■ V)u (3.29) 

4. Energy Principles 

We shall now show that the energy of steady flows is stationary compared to the energy 
of any nearby trial configuration with the same total mass, linear and angular momentum 
and the same mass and vortex flux in displaced fluid elements. 

4-1- Some Differential Identities 

Let F be a function of {a,P,iJ,) = a''. Since 

= = • ^^^^= (4.1) 

it follows that 

A(W) = -V^^dka^ = -VC • Va' (4.2) 
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and therefore 

AVF = -V^^^feF + VAF = -Vf • VF + VAF (4.3) 
In particular, see (3.9), 

A^^ = aAVP + AVu = -Vi- W + VAu 

= VW -i+ViAu-W -i) (4.4) 

Where Au is defined by varying (3.10) 

A[V • (pW)] = AV • [p{-aVP + ffc)] (4.5) 

This defines Au by an eUiptic equation of the same form as v itself 

V • [pV(Azy)] = etc... (4.6) 

and therefore Au is defined with the same type of Green function as u. Another needed 
variation is that of $ which, according to (2.1) and (3.5) may be written 

^ = -g[ = -Gf'^ (4.7) 

"^J \f{a)-f'{a')\ "'i R ^^■'> 

Thus 

A^ = -G j [i- v| + i- V'^ d'a' (4.8) 
and therefore the variation of the gravitational potential energy 

A J ^^pd^x = aJ ^^d^a = J ^A^d^a = j ^^pd^r. (4.9) 
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4-2. The Variation of the Energy 
The energy of the flow 



E 



-I 



p(fx (4.10) 



in which e{p) is the specific internal energy of the barotropic fluid, related to the pressure 
and the specific enthalpy: 

e{p) = h-- P = -P'^ (4.11) 
p dp 

With AW given in (4.4) and with A J ^^p(fix in (4.9), one readily finds the following 
identity for AE, 

^E = J {[{W ■ W)W + V(/i + $)] • |p - V • {pW){Aiy - W ■ i)}(fx 

+ y V • [{Au - W ■ i)pW - Pi](fx. (4.12) 

4-3. The Constrained Variational Identity and the Principle of Stationary Energy 
We can now also calculate AP and A J which according to (2.11) are given by 

AP = J AWd^a and AJ = J ^ x W + f x AW)d^a (4.13) 
Using (2.3) that defines U, we then obtain a new identity: 

AE -b- AP -n^- AJ = J [d-ip + U{W ■i-Ai^)]d^x 

+ [ [{Aiy - W ■ i)pW - Pi] ■ dS, (4.14) 
Js 

where O and U have been defined in (2.2) and (2.7), respectively. 
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We see from (4.14) that if p\s = P\s = 0, mass is preserved {U = 0) and Euler's 
equations for steady flows in uniformly moving coordinates hold {O = 0) then the energy 
is stationary {AE = 0) when linear and angular momenta are kept fixed (AP = AJ = 0). 

4-4- Weak and Strong Principles of Minimum Energy 

Reciprocally if AE = with AP = A J = (which defines Ab and Al^c), then the 
following must hold: 

(a) Either v is defined by mass conservation and the principle of stationary energy 
provides Euler's equations = 0, 

(/?) Or u is an additional variable and AE = if, in addition to O = 0, we have mass 
conservation, U — 0. 

The principle of energy with four instead of three independent functions, r and v, 
we shall call the weak energy principle, as opposed to that with only three independent 
functions which we call the strong energy principle. 

As far as stationarity of energy is concerned, the distinction is of no consequence. 
But it becomes important when we consider second order variations and the principle of 
minimum energy. Second order variations of (4.14) near the extremum AE = are given 

by 

A'^E - b- A^P - flc ■ aV= J [AO ■ Cp - A2A{Av - W ■ i)]d^x (4.15) 

Notice that with (2.9) and (2.10), the boundary term in (4.14) does not contribute to A'^E. 
The relations A^P = A^ J = define A% and A'^Vlc. Stationary configurations are stable 
if 

A'^E > 0. (4.16) 
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The strong energy principle reads then 

^''E strong = J ' ^P*^^"^ > ^ (^-l^) 

which contains Au that must be obtained from (4.6). The weak energy principle is 

^^E^eak = j[^d- ip - AU{Au - W ■ 0]d^x > (4.18) 

and is manifestly simpler to apply since Av is here independent. The weakness of this 
principle compared to (4.17) can be understood by the fact that variations of the trial 

— * 

functions are mass preserving and isocirculational but that the velocity trial field W does 
not conserve mass. Only extremal W's do. Thus fluctuations in A'^E would include 

— * — * — * 

Al^^'s that do not necessarily satisfy the mass conservation equation. W + AW being less 
restricted than a real dynamical + Al^, some instabilities might show up that cannot 
exist and would not show up in a perturbation analysis. For this reason this energy principle 
is weaker. But the weak principle may be helpful when the strong principle involves too 
hard calculations. Moreover, in numerical calculations, in which the computer "searches" 
the minimum of E, one has to find the Green functions G{f, r') for v at every step and that 
may be time consuming. Therefore the energy principle with four independent functions 
may indeed be useful. 

5. Principles of Stationary Action 

5.1. Introduction 

It may appear to make little sense to set up a Lagrangian formulation in terms of 
variables that take a priori account of all conservation laws of motion. First the dynami- 
cal equations insure automatically that the constants of the motion are indeed constant. 
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Second, in hydrodynamics, fixation of the values of constants of motion does not reduce 
very much the number of independent variables but complicate considerably the equations 
of motion. 

However, it is important to show that our labelling defines proper Lagrange variables, 
even if we are never going to use them, for the following reason. An ordinary Lagrangian 
has the form L = T — V while the energy E = T + V {T the kinetic energy, V the 
potential energy). If T is purely quadratic in the time derivatives of the Lagrange variables, 
A^E > is a necessary and sufficient condition of stability. It is however well known that 
Lagrangians subject to non-holonomic constraints, like circulation conservation or fixed 
linear momentum, contain additional terms, so called gyroscopic terms G , that are linear 
in time derivatives. The form of L is then T + G — V, the energy, however, is still T + V. 
In those circumstances, A'^E > is only a sufficient condition of stability, not a necessary 

— * — * 

one. In our representation in which not only fixed J and P but also fixed circulation and 
mass conservation have been incorporated, we may have lots of gyroscopic terms in the 
Lagrangian. It is therefore useful to consult the form of L and to see if there are gyroscopic 
terms so as to know when our energy principles give necessary and sufficient conditions of 
stability or only sufficient ones. For this reason we have first to proof that the Lagrangians 
of our flows provide indeed Euler's equations. 

5.2. Time Dependent Lagrange Variables and Dynamical Equations 

We use the same Lagrange variables a'^ introduced in section 3. Now, however, they 
depend also on the time: a'' = a''{f,t). Surfaces of constant a'^ define a velocity w which 
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satisfy the equations: 

a''{f+wdt,t + dt) = a''{f,t); (5.1) 

in the hmit dt ^ 0: 

a'' +w-Wa'' = 0. (5.2) 

Since arbitrary displacements of a fluid element with labels (see section 3.4) are both 
mass preserving and isocirculational it follows from (5.2), (3.4) and (3.6) that mass and 
vortex strength are both preserved along a motion with velocity w: 

p + V • (pw) = (5.3) 

and 

uj + rot{Co X to) = (5.4) 

There exists a useful explicit expression for w that simplifies various formulas; consider 
the following identity obtained from r(t, a^): 

r = r[t,a^{r,t)] (5.5) 

Since the right hand side must be independent of t, we have: 

extracting from (5.2) and replacing it in (5.6) gives then: 

Br 

{ttX = w (5.7) 



dt 

and, quite generally. 



{^^^)a = F + w-WF (5.8) 
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The "absolute" velocity of the fluid is, say, W plus the velocity of the vortex lines w: 

iJ=^ + W (5.9) 

The relative velocity is similarly: 

u = w + U (5.10) 

W and U are related by eq. (2.3). The equation of mass conservation along the motion of 
the fluid must be 

p + V-{pu) = (5.11) 
and (5.3), with (5.10) and (5.11), imply that U satisfles equation (2.7) again: 

U = V ■ (pU) = 0. (5.12) 

The dynamical equations of motion are slightly different from (2.2) and are well known: 

Od {u-V)v + nc X w + V(/t + $) = 0, (5.13) 

This familiar equation may also be written in a slightly less familiar form which will be 
useful; with (5.8) and (5.10), we have 

Od = {-^)a + {U ■ V)v + VL^ X ^7+ V(/i + $) = 0, (5.14) 

In stationary flows w = Q ,v = W and (^^)a = ; Od becomes then the O of (2.2). 

5.3. The Action 



The Action of the system is given (see KLB 85) by: 

'to 
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j\dt = J J[^w^ - {^W^ + e+^^)]pd^xdt^ J J[^w^pd^x-E{t)]dt (5.15) 



Notice the minus sign and the fact that E appears here effectively as the potential. The 
Lagrangian of isocirculational flows is actually a Routhian that incorporates mass and 
circulation conservation. With (5.9) we may write L as this: 

L = J[w -v- (^^n^ + e+^<^^]p(fx (5.16) 

which is of great help in calculations. The Lagrange variables are thus {a, (3, /x) = as 

— * 

in stationary flows; v, instead of W, has now a Clebsch form: 

v = aVP + Vi^. (5.17) 

where is a single valued functional of deflned by (5.12) and is linear in {a^,b,(lc) 
because w is linear and homogeneuos in a'^. Quantities like e and $ depend on p only 
which is a functional of a''. The Action is thus a functional of , quadratic in but 
not quadratic homogeneous. Notice that L is apparently linear only in through w. 
However, L is quadratic in u which is also linear in w, b and Qc- The constants b and fie 
are deflned by the values of linear and angular momentum which are now given by: 

P = Jvp(fx = J{aVP + Vu)pd^x = 0, (5.18a) 

and 

J = Jrxvpd^x = jrx{aVI3 + Vv)pd^x = jQ\ (5.186) 

— * — * , 

b and Q,c are thus equally linear in a through v. 

5.4- Variational Identities for Calculating A A 

Let us now make small displacements ^ of fluid elements and calculate the changes 
A^, AP and A J in A, P and J . Much of the calculation has actually been done in 
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section 4. Indeed notice first that E{t),P and J are simifiar to E in (4.10), P in (2.11b) 
and J in (2.11c) with v replacing W. Since v, in time dependent flows, and W, in steady 

— * — * — * 

flows, are both represented by aV/J+Vi^ , we obtain, straightaway AE{t) — b- AP — Qc- 
U by u: 

AE{t) - b- AP -Qc- AJ 

= jm . v).+ X + n • 6 + V . (p.)(. . r- a.)}.3. 

+ / [{Av - V ■ Opu - Pi\ ■ dS, (5.19) 
Js 

Notice that (5.19) is a variational identity; v has been treated as an independant function. 

To calculate there remains, however, some work to do on the v ■ w part of L. Following 
(5.7) and (5.17) 

A{w-v) ^ A{—)a ■ V + w ■ A{aV (3 + (5.20) 
We can use (4.4) to rewrite A{a'V(3 + Vz^)and obtain 

A{w-v) = {-^)a ■ V + {w ■ V)v ■ ^ + w ■ V{Au - V ■ ^) (5.21) 

so that, with some integration by part and, with the help of (5.8), we can write 

A J w-vpd^x= ^(/^^^^^^^) ^ j [_i/.fp+V-(pt(;)(i;-|'-Azy)](i3a; + j {Aiy-v-i)pw-dS, 

(5.22) 

Thus, J[AL + b- AP + Qc ■ AJ\dt is given by (5.22) minus (5.19); in writing the result we 
shall use Od deflned in (5.13) and U in (5.12): 

A^ + J{b-AP + nc-AJ)dt={Jv-ipd^x)\ll-J j[dD-ip + U{v-i- Av)]d^xdt 
[{Av - V ■ i)pW - Pi] ■ dSdt, (5.23) 
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5.5 Weak and Strong Principles of Stationary Action 

Identity (5.23) leads straightaway to the following results. If p|s = P\s = 0, mass 
is preserved (W = 0) and Euler's equations hold {Od = 0) then the Action is stationary 
(AA = 0) when linear and angular moment Reciprocally, if AA = with AP = A J = 
(which define A6 and Ailc ) then the following must hold: 

(ck) Either v is defined hyU = Q and the principle of stationary Action provides Euler's 

— * 

equation = . 

{(3) Or v is an additional variable and A is stationary if in addition to Od = 0, we 
have mass conservation U — 0. The principle of stationary Action, with four independent 
functions, and instead of three, we call the weak principle of stationary Action. It is 
the dynamical counterpart to the weak energy principle. 

6. Gyroscopic Terms for Small Deviations from Steady Flows 

6.1 The Action for Small Deviations 

Small time dependent deviations from a stable flow are given by Euler's linearized 
equations whose Action is the second order variation A'^A of A calculated at the stationary 
"point" . Let fo{t, a'^) be the coordinates of the fluid element of a particular stationary flow 
and fo + |(t, a'') the coordinates in the perturbed flow. The gyroscopic terms in A'^A are 
the bilinear functionals of ^ and ^. The gyroscopic terms of the fully constrained Action 
are given by the I*, ^ terms in J[A^L + 6- A'^P + Qc ■ A'^J]dt as evaluated at the stationary 
point. Let the 

Wo = 0, vo = Wo, Uo = Uo, V • (poUo) = etc... (6.1) 
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Notice that: 

(§)a|o = f+(^o-V)f=e' (6.2) 

With (6.1) and (6.2) we readily obtain A^L + etc... from (5.23) at the stationary "point". 
Notice that because of (2.9) and (2.10), the second variations have no surface terms. Since 
the terms at t = to or ti do not contribute either, we have: 

[A'^A + J^\b- A^P + Qc- ^^j)dt]o = J J[-AdD-ip + AU{Ai^-v-i)]odtd^x (6.3) 

in which 

AU\o = A[V • (pf7)]|o = V • [po(VAzy - (vf • Uo + Uq- Vf) - • ff 

where ffc is the vector defined in (2.3); of AOd we write only those parts susceptible to 
contribute to the Gyroscopic term: 

( J J -AOd ■ ipodtd^x)Gyro = J J {(■ {VAu - • Wq) + [(f - VAu + A-XXX 

6.2 Gyroscopic Term for the Strong Principle. 

In the strong energy and Action principle, v is defined by W = and Av is defined 
by AU\o = 0. Av may be decomposed into a "dynamical contribution" Avd which is zero 

— * 

when ^ = and a "steady" part Avg: 

Av = Avd + Avs (6.6) 

Similarly, ry^ has a dynamical and steady contribution (see our remarks at the end of 
section 5.3) and we write 
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Arfc = AficD + Affcs (6.7) 
AU\o = gives the following two equations for Avd and Ai^s deduced from (6.4) 

V • (poVAud) = V • [po(e + AffcD)] (6.8a) 

V • {po^Aus) = V • [po{Uo ■ Vf + Vi- Uo) + poVe • r/c + Affas)] (6.86) 

Gyroscopic terms in (6.3) will only come from AOd since AU\o — 0. Using eq. (6.8) 
the bilinear terms of [A^^ + J^^ib- A^P + Qc ■ A'^J)dt]o, which we denote by A'^Gstrong 
are 




A^G strong = J J [^Aud ' VAug - C " • Wojpod-'xdt 

[AffcD ■ ^Aus + Affcs ■ VAvD + ArycD • V^- #o 
+ {AVt^D X #0 + X VAvd) ■ ilpod^xdt, (6.9) 

This expression becomes somewhat simpler for motions that are steady in inertial 
coordinates {ffc = 0), when, for some reason, Az/c does not contribute either: 

A^Gstrong ^ J J i'^^^'D ' ^Aus - f Wo]pod^xdt (6.10) 

We conclude that the strong energy principle provides a necessary and sufficient condition 
of stability if: 

A'^Estrong > and A^Gstrong = (6.11) 
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6.3 Gyroscopic Term and the Weak Energy Principle. 

The weak Lagrangian is linear in time derivatives; L is of the form G — V since T = 0. 
The perturbed Action has no quadratic term in ^ which appears only in /S^Gyj^ak- As a 
result, if A^G^^eafe 7^ 0, A^iJ^gafe > is certainly a sufficient condition of stability (this 
can easily be shown). However, if /Sp'Gy^eak = and /Sp'Ey^eak > 0, the linearized equation 
of motion have no solution at all! Usual arguments about stabili 

It is than still interesting to obtain A^G^^eafc ^ot., however, very different reasons than 
we wanted IHi^ G strong- In the weak Action principle, P and J do not define h and Q,c 
[see (5.18)] because v is independent of them. Instead of describing the motion in moving 

— * — * — * 

coordinates, we rather stay in inertial coordinates and fixP=J — Jo = with Lagrange 
multipliers, say, h{t) and VLc{t). The Action is then 

= A + J [b{t) ■ P + ndt) ■ {J- Jo)]dt (6.12) 

Among all the dynamical perturbations, we consider only those for which P = J — Jq = 0. 
The formal expression for A^A^ is than exactly the same as the right hand side of (6.3). 
The gyroscopic term A^G^eak 

(A2G)„eafc = J J [(21'+ ^Vcd) ■ (VAu - Ve • Wo) - A^cD X Wo ■ ^pod^xdt (6.13) 
We conclude that a sufficient condition for stability is 

A^E^eak > and A^G^eak ^ (6.14) 



This work was started when S.I. visited the Racah Institute of Physics, Hebrew Uni- 
versity of Jerusalem in 1987. He acknowledges its support and hospitality. 
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FIGURE CAPTIONS 

;. 1. Representation of a trial configuration indicating a surface of constant load A = Aq, 
a cut of constant (3 = (3q hanging on the central "string" and the line A = 0. The 
point P, at a shortest distance from the 2;-axis helps to define the orientation of the 
x, y plane. 
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